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ABSTRACT 

The most accurate way to measure the energy levels for the O II 2p 3 ground 
configuration has been from the forbidden lines in planetary nebulae. We present 
an analysis of modern planetary nebula data that nicely constrain the splitting 
within the 2 D term and the separation of this term from the ground 4 S3/2 level. 
We extend this method to H II regions using high-resolution spectroscopy of the 
Orion nebula, covering all six visible transitions within the ground configuration. 
These data confirm the splitting of the 2 D term while additionally constraining 
the splitting of the 2 P term. The energies of the 2 P and 2 D terms relative to the 
ground ( 4 S) term are constrained by requiring that all six lines give the same 
radial velocity, consistent with independent limits placed on the motion of the 
+ gas and the planetary nebula data. 

Subject headings: atomic data — methods: data analysis — techniques: spectro- 
scopic — ISM: individual (Orion Nebula) — planetary nebulae: general 



1. Introduction 



In this paper we determine the four energies which describe the separation of the five 
ground configuration energy levels of O II shown in Figure 1, through analysis of observa- 
tions of the visible forbidden transitions. Being forbidden transitions, these are weak and 
not accessible in the laboratory. The ground configuration levels can be determined from 
laboratory data on permitted ultraviolet transitions from higher states, but in practice the 
uncertainties in the derived energy levels are larger than those obtained from the nebular 
visible transitions. Thus the published compilations of data on these energy levels (Eriksson 
1988; Martin, Kaufman & Musgrove 1993) rely heavily on the pioneering work of Bowen 
(1960) and de Robertis, Osterbrock, & McKee (1985) on planetary nebulae (see Table 1). 
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Our analysis of new data improves upon these energies and provides a discussion of the 
uncertainties. 

Bowen (1960) observed the two blue [O il] lines in seven bright planetary nebulae (a total 
of sixteen photographic plates) and the four red [O n] lines (blended into two lines) in one 
nebula (three separate plates). The + gas of the planetary nebulae was assumed to have an 
average bulk velocity identical to the measured velocity of the H + gas (the H + energy levels 
being well known). The observed [O il] wavelengths adjusted to this velocity give absolute 
wavenumbers for the energy level differences for the transitions. These early data resolved 
the 2 D term into its two J components, but were not able to determine the splitting of the 2 P 
term. The latter splitting was determined in de Robertis et al. (1985) using high- resolution 
digital spectra of the four red lines in one bright planetary nebula (NGC 7027). As can be 
appreciated from Figure 1, the red transitions ending on a common lower level constrain 
the splitting of the 2 P term while the transitions from a common upper level constrain the 
splitting of the 2 D term, independently of the blue data, de Robertis et al. (1985) found a 
value for the latter splitting somewhat different than that found by Bowen. Since absolute 
velocities or wavenumbers were not measured, de Robertis et al. (1985) calculated the 2 D- 
2 P separation by assuming that the A7319/A7320 line blend measured by Bowen (1960) was 
actually a measure of the stronger A7320 line, and that Bowen (1960) 's roughly symmetric 
A7330/A7331 blend was a measure of this line pair's average. 

We show that H II regions can also be used to determine the energy levels of O II. With 
high-resolution echelle spectroscopy of the Orion nebula (resolving all red and blue lines), 
we are able to measure not only the splitting of both the 2 P and 2 D terms but also the 
separation between these two split terms. The observations are described in § 3 and the full 
analysis, including uncertainties, in § 4. But first we update Bowen's work by analysis of 
more recent digital data on 23 planetary nebulae (§ 2). This provides important insight into 
the method of analysis of the full Orion data and also provides improved values of the two 
energies constrained by the blue lines a la Bowen. 

We conclude by comparing our results with recent work by Sharpee et al. (2004). 

2. Bowen revisited 

Using a coude spectrograph on the Hale telescope, Bowen (1955) observed the blue 
spectra of seven planetary nebulae photographically. Since the introduction of the Hamilton 
Echelle Spectrograph in 1987 (Vogt 1987), high-resolution spectra of (at least) 24 planetary 
nebulae spectra have been published, the relevant data for our purposes being the line lists 
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with accurately tabulated wavelengths. (One nebula, NGC 6818, had to be discarded because 
of a grossly discrepant wavelength - we suspect a typographical error.) These new spectra 
resolve the blue lines but not the red pairs of lines. 

Following Bowen we therefore have two lines to constrain two energies. But there is 
an important requirement: we need to know the velocity of the line-producing gas and any 
uncertainty (Svo+) will have predictable consequences. For example, if one uses the 3729 A 
line to find the energy difference for this transition, -E" 4 s 3 / 2 - 2 d 5/2 (note the convention of an 
"upwards" transition), because of this degeneracy the Doppler effect gives an uncertainty 

<^s 3/2 -3d 5/2 = 0.0895^ o+ , (1) 

where throughout this paper energies are in units of cm -1 and velocities in km s _1 . Looking 
ahead to the red data, 

^D 3/2 - 2 p 3/2 = 0.0455^o+ = 0.50865^4 S3/2 _ 2D5/2 , (2) 

but note that the smaller splittings i& D _2 D and £ , 2p 3/2 _2 Pi/2 are much less susceptible 
to any uncertainty in the velocity. 

In Bowen's work, permitted lines of H + and He + gas set the rest-frame velocity and 
similarly here all wavelengths for all 23 planetary nebulae are first put into such a reference 
frame with H + at rest. As tabulated by the authors, the planetary nebulae wavelengths 
have in fact already been corrected for the previously-known systemic velocity of the whole 
nebula, but nevertheless we find from the data that the H + is apparently not quite at rest. 
This is either because the systemic velocity used was only approximate, or because for the 
slit positions observed the mean velocity of the H + gas is not identical to that averaged 
over the nebula, which would not be surprising given incomplete coverage of an expanding 
nebula. To define the H + frame for each nebula, we used the eight or nine unblended H I 
lines (H16 to H<5) near the blue [O il] lines, these all being contained in the same echelle 
spectrum. The corresponding blue wavelengths are given in Table 2. 

The measured velocity of + is not necessarily the same as that of H + because of the 
ionization structure and expansion that exist in the nebula and the fact that the slit does 
not usually cover the entire nebula. Bowen's approach was to average results over several 
nebulae, since on average the two velocities should be equal. Given data on 23 nebulae, we 
can improve upon this iteratively as explained below. 

Our analysis is basically to develop a parameterized model, and then optimize the 
parameters by non-linear least squares to match predicted wavelengths, in air, with the 
wavelengths tabulated. The parameters of the model are the energies and any velocity 
offsets (vq+) deemed necessary. (Toward this end, deviations of observed wavelengths from 



-4- 



the model predictions are expressed in terms of velocity.) The energies can be taken as 
the successive energy differences, four independent values in the full model, or simply the 
energies of the four upper levels. Even with the first of these two options there is covariance 
in the resulting solution, since four of the six energy transitions (corresponding to four of the 
six available wavelengths) couple the independent successive energy differences. The "model 
uncertainties" from the goodness of fit to the model are the 68.3% confidence intervals for 
one-dimensional marginal distributions for each of the parameters. 

Let us return to the blue lines, for which we have 46 measured wavelengths for 23 
planetary nebulae. As energy parameters, we used E4 Ss 2- D& an d E 2 d 5/2 ~ 2 d 3/2 - m the 
initial model we went to the extreme of introducing 23 velocity offsets. This precludes 
determining £74s 3/2 _2d 5/2 and we find £ , 2 D5/2 _2 D3/2 = 19.79 ± 0.05. For each nebula the two 
velocity residuals (from the two wavelength residuals) are of equal magnitude (denoted crj nt ) 
with opposite signs, indicating that relative to the model the two lines are too close together 
or too widely separated. Overall the rms velocity residual was 0.96 km s -1 . Those nebulae 
with considerably larger rms values can be judged to have data of lower quality; i.e., even 
with the luxury of the maximal number of parameters, the data are still not going to be 
well matched by the model. Two nebulae with residuals greater than 2.4 km s _1 have been 
assigned weight 0.25 (in the calculation of x 2 ) while another four with residuals greater than 
1.2 km s -1 have been assigned weight 0.5. No bias is introduced in subsequent calculations 
of the splitting -E 2 d 5/2 - 2 d 3/2 since equal numbers of "too close" and "too separated" cases 
are involved; now £ , 2 D5/2 _2 D3/2 = 19.79 ±0.04. The weighted rms residual is 0.71 km s" 1 and 
no subsequent model, with fewer parameters, can improve upon this. 

The next model goes to the other extreme, fitting only the two energy differences, finding 
£4s 3/2 _2 D5/2 = 26810.68 ± 0.13 and £ 2D5/2 _2 D3/2 = 19.79 ± 0.18. The latter energy difference 
is still close to that in the initial model, but is determined with less confidence because the 
two-parameter model fits the data less well. The rms velocity residual is 6.7 km s _1 , much 
larger than suggested by our assessment of the data quality, and so clearly indicating a less 
than optimal model. 

Thus our goal was to improve the model iteratively by adding a minimal number of 
parameters Vq+, for a subset of the nebulae, expecting a significant reduction in y 2 per 
degree of freedom. We identified seven nebulae for which the residuals exceeded 6.7 km s _1 
and also 3<7i n t < 6.7 km s _1 and for these included a velocity offset parameter. This produced 
a markedly improved model with an rms residual of 1.9 km s _1 and -^ 4 s 3/2 - 2 d 5/2 — 26810. 77± 
0.04 and -E ,2 d 5/2 - 2 d 3/2 — 19.79 ±0.06. Repeating this process with the new rms identifies eight 
more nebulae which would benefit from velocity offsets, for a total of 15 of the 23 nebulae. 
This iteration reduces the rms velocity residual to 0.82 km s" 1 ; this is now comparable to the 
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above estimate of the quality of the data, indicating that adding further parameters would 
not be justified. £ 4S3/2 _2 D5/2 = 26810.77 ± 0.03 and £ 2Db/2 _2 D3/2 = 19.79 ± 0.03 - the only 
difference with respect to the previous iteration being a lowering of the error (confidence 
interval). It is important to acknowledge that there is a systematic error in -E ,4 s 3/2 - 2 d 5/2 
of order 0.03 cm -1 because the introduction of parameters t>o+, while hopefully unbiased, 
is still subjective. Such systematic errors are recorded separately in Table 1 to distinguish 
them from confidence intervals. 

The energies determined are close to those given by Bowen, but significantly different 
than those derived by de Robertis et al. (1985) - refer to Table 1. 

3. Observations of the Orion Nebula 

High-resolution echelle spectra were obtained over the course of two nights in 1997 (5100- 
7485 A) and 1998 (3510-5940 A) with the CTIO 4-m Blanco telescope (refer to Baldwin 
et al. (2000) for details). Three different lines of sight were observed, referred to as 1SW, 
x2 (Baldwin et al. 1996; Rubin et al. 1997) and 37W (Baldwin et al. 1991). To extract the 
information on individual lines, in particular the wavelength, the data were modeled with a 
Gaussian (three parameters: central wavelength, FWHM, area) and a linear baseline (two 
further parameters). The rest-frame velocity of the H + gas along each line of sight was 
determined from the six strongest unblended H I Balmer lines. 

All six [O n] forbidden lines are seen with good signal to noise in these spectra. Data 
for 1SW are shown in Figure 2. The pairs of [O il] red lines were slightly blended and so 
were analysed using a double Gaussian fit. The results of our line fitting are summarized in 
Table 3. All line profiles are similar as seen in the matching FWHM and directly from the 
spectra in Figure 2. 

For the common upper level line pairs A7320/A7331 and A7319/A7330 the line strength 
ratios can be predicted directly from the transition probabilities (Zeippen 1987; Wiese et al. 
1996), offering an independent check of one aspect of the fits. Results are presented in 
Table 4. There is reasonable agreement between the theory and the observations. The sole 
anomaly, seen in the A7319/A7330 ratio in the x2 line-of-sight, arises because of a velocity- 
shifted component from a photoionized Herbig-Haro shock (Blagrave & Martin 2004a) that 
has 2-4% of the nebular [O il] flux. The nebular A7330 and A7319 lines are contaminated 
by the velocity-shifted components of A7331 and A7320, respectively. A higher relative 
contamination from the stronger of these two lines, A7320, results in a higher A7319/A7330 
ratio. 
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All four [O n] red lines are found in the same echelle order of a single exposure, unlike in 
the data of de Robertis et al. (1985) where the A7319/A7320 and A7330/A7331 line pairs were 
obtained in two separate spectra. Measuring the 2 P energy splitting £- 2 p 3/2 - 2 p 1/2 depends 
on the wavelength difference within each of the line pairs and so the results in de Robertis 
et al. (1985) should be accurate. On the other hand, the splitting of the 2 D term fi D 2 D 
depends on the wavelength difference between the pairs and so our single spectrum containing 
both line pairs, with only a single wavelength calibration in the same echelle order, should 
yield a more accurate result. 

All lines produced by the same ionized species should have the same velocity. However, 
using the current published energy levels and derived rest wavelengths in air (Eriksson 1987; 
Martin et al. 1993) together with our observed wavelengths, we obtain + velocities that 
are grossly inconsistent, well beyond the uncertainties propagated from the measurement 
errors of the observed wavelengths (see Table 5). In particular, there is no explanation why 
lines in the same wavelength region and originating from a common upper level ( 2 Pi/2 or 
2 ?3/2) should yield significantly different velocities, as is observed to be the case in columns 
7 and 8 of Table 5; at the very least, there is a problem with the splitting -E ,2 d 5/2 - 2 d 3/2 - 
The consistency of the data for the three lines of sights, and the lack of agreement of the 
velocities from all of the lines, points to inaccurate [O n] rest wavelengths arising from 
poorly determined ground configuration energy levels, including the separations i?4 S _2 D 
and i?2 D _2 P . This was the original motivation for this paper. 



4. Constraining the Energy Levels with the Orion Nebula Data 

With a set of six accurate wavelengths, for each of three lines-of-sight (1SW, x2, and 
37W), it is possible to obtain the energies of the four excited levels in the ground configuration 
of O II (Fig. 1); the problem is over-constrained. However, as encountered in the analysis 
of planetary nebula spectra in § 2, there is the possibility of an unknown velocity offset of 
the + gas for each position observed. But even with an extra velocity offset parameter v + 
for each position, the problem is still over-constrained (even for a single position) and thus 
amenable to modeling and least-squares optimization. 

As with the planetary nebulae, there is a Doppler-related degeneracy to be resolved as 
well, through independent constraints on v Q + . We show how this is possible in Orion, it not 
being sufficient to assume that on average, over many positions, the velocities of the + gas 
and H + gas are identical (in which case (v +) = 0). 

The Orion nebula can be represented by a blister model where gas is accelerating toward 
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the ionizing star and the observer away from the background molecular cloud (Balick et al. 
1974). Because of ionization stratification in this accelerating flow, the velocity is correlated 
with the ionization potential (I.P.) : gas that is more highly ionized is more blue shifted (has 
more negative velocity). This is observed; see, e.g., Baldwin et al. (2000). Figure 3, plotting 
velocity against the emitting species' LP., summarizes this effect for our observations of 
many lines of many different ions for the three lines of sight. Models of the nebula (Baldwin 
et al. 2000) show that the + zone is relatively narrow (see their Figure 4) and so there 
should be a well defined velocity. From the trend seen in Figure 3, the expected value for 
v +, at l.P. 35 eV for + , lies between and 5 km s _1 , set in part by the [S ill] (l.P. = 
34.79 eV) velocities at similar ionization potential. Plotting the results in columns 7 and 8 
of Table 5 in Figure 3 would clearly reveal their discrepancy, again pointing to a problem 
with the energy levels. The velocities from our new + model clearly satisfy the general 
constraint set by the surrounding lines. 

As in the analysis of the planetary nebula data, our examination of the Orion nebula 
data is rooted in a model of the energy levels, used to predict the air wavelengths. In this 
case we set £' 4 S3 /2 - 2 d 5/2 — 26810.77, the result already found in § 2. The remaining analysis 
is then a non-linear least squares fit (unweighted), using six parameters in the model: the 
other three energy differences, and the three velocities vo+ ■ This model produces rms energy- 
difference residuals of only 0.4 km s" 1 . Line by line and position by position, this model 
gives the velocities listed in column 9 of Table 5; these are clearly now quite consistent with 
one another. 

This choice of -E ,4 s 3/2 - 2 d 5/2 produces vo+ equal to 4.0 ± 0.3 km s _1 , 4.1 ± 0.3 km s _1 , 
and 1.4 ± 0.3 km s" 1 for 1SW, x2, and 37W, respectively (because of changes in geometry 
from one line-of-sight to the next, these velocities need not be identical). It turns out that 
these fall in line with the trend in Figure 3 and so no optimization was carried out on this 
energy difference. 

Recall however, that if a slight change were made in this energy, t> + would respond ac- 
cording to equation 1. Thus the model uncertainty <5_E4g 3/2 _2 D5/2 from the planetary nebulae 
analysis, ±0.03 cm -1 , results in Sv + = 0.34 km s -1 . This is of the same order as the model 
fitting uncertainties for the Orion data. In principle, one might start by determining Vq+ 
and its uncertainty from interpolation in Figure 3 and work backwards to E4 S _2 D and 
its uncertainty. Given the consistency, we adopt the tighter constraint found independently 
in § 2. This will propagate as a systematic uncertainty through to the other energies (see 
equation 2). 

Uncertainties in vq+ have the same effect as uncertainties in the wavelength calibration, 
which we deduce is accurate to 0.5 km/s, consistent with the scatter of individual blue 
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H + line velocities about the mean. Thus, for the remainder of the Orion analysis, we adopt 
Svo+ = 0.5 km s _1 , noting that this subsumes errors in the wavelength calibration of the blue 
lines. For the red lines there is an additional systematic error of order 0.15 km s _1 relating to 
the alignment of the separate red and blue echelle spectra into the same velocity /wavelength 
system. The next issue is the accuracy of the wavelength calibration for individual lines. 
To assess this, we note that there are numerous lines that are duplicated in neighbouring 
echelle orders of the same spectrum. By comparing the duplicate measurements of these 
lines' wavelengths, we find a residual characteristic order-to-order difference of 0.65 km s _1 , 
independent of measurement uncertainty. In practice, we feel that measurements of the red 
[O n] lines, which all appear in the centre of an order, have a systematic uncertainty of less 
than half this, ~ 0.3 km s -1 . Taken all together, the combined systematic error is of order 
0.7 km s _1 , which corresponds to 5E2 D _2 P = 0.03. Therefore, the separation -E 2 d 3/2 - 2 p 3/2 
is properly quoted as 13637.34 ± 0.01 ± 0.03. This is significantly different than what has 
been previously adopted (see Table 1). 

On the other hand, for the small splittings, the systematic effects are tiny compared to 
the model uncertainties, and so are listed below as 0.00. As anticipated, we find that the 
value of the splitting Et-p 2 P of the 2 P term, 2.02 ± 0.01 ± 0.00, is in close agreement 

3/2 1/2 

with that found by de Robertis et al. (1985), 2.00 ±0.03 (largely adopted by Eriksson (1987) 
and Martin et al. (1993)). 

The splitting -E ,2 d 5/2 - 2 d 3/2 of the 2 D term is 19.80 ± 0.01 ± 0.00, closely consistent with 
the value 19.79 ± 0.03 ± 0.03 obtained in § 2 from only the blue lines of planetary nebulae. 
We conclude that the splitting of the 2 D term should be revised from the value 20.1 ± 0.1 
found by de Robertis et al. (1985) and (largely) adopted by Eriksson (1987) and Martin 
et al. (1993). 

These energy differences allow us to calculate the energies of all four energy levels (see 
Table 1). The model was rewritten with these energies as the parameters in order to track the 
effects of covariant changes in this non-linear model of the air wavelengths, and thus provide 
the appropriate one-dimensional marginal confidence intervals reported in Table 1. The 
errors reported are again the model fitting uncertainty and the systematic error. The latter 
reflects the propagation of the error in -E ,4 s 3/2 - 2 d 5/2 phis the systematic error in -E' 2 d 3/2 - 2 p 3/2 
for the highest two energies 

These energy levels allow for the calculation of the air wavelengths of lines in the UV in 
addition to those in Table 5. From £4s 3/2 - 2 p 3/2 E* s 2 P we calculate A2470.347± 0.001 ± 
0.003 and A2470.223 ± 0.001 ± 0.003, respectively. 

Because of the presence of systematic errors, the above analysis was carried out with 
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equal weights for all lines. However, we have repeated the analysis using a weighted fit based 
on the uncertainties found in measuring each line. -E 2 d 5/2 - 2 d 3/2 increased by only 0.001, 
-^ 2D 3/2- 2 P3/2 decreased by 0.015 (well within the systematic uncertainty) and ^ 2 p 3/2 - 2 p 1/2 
increased by 0.007 (recorded in Table 1 as 0.01 systematic uncertainty). 

5. Comparison to Recent Work of Sharpee et al. 

Sharpee et al. (2004) have observed the four red lines (AA7319, 7320, 7330, 7331 A) 
in sky spectra using the high-resolution echelle spectrograph (HIRES) on Keck I. In their 
determination of the O II 2p 3 energy levels, they also make use of data-sets from nebulae 
(Bowen 1960; de Robertis et al. 1985; Baldwin et al. 2000; Sharpee et al. 2003). The nebular 
data are given lesser weights (0.018, 0.16, 1.00, 1.00, respectively) than the HIRES data-set 
(10.00). Velocity corrections for data-sets of three nebulae (de Robertis et al. 1985; Baldwin 
et al. 2000; Sharpee et al. 2003) were found using a weighted least-squares analysis similar 
to what was done in § 2; Bowen (1960) is an average of multiple nebulae and was assumed 
to have (v +) = 0. From these five data-sets, they obtain splittings of 19.810 ± 0.006 and 
2.010 ± 0.005 for the 2 D and 2 P terms, respectively - both in excellent agreement with our 
results. 

The 2 D 3/ 2- 2 P 3 /2 separation, 13637.403 ± 0.004, is quoted here with their la model 
fitting uncertainty. Even though for the HIRES data one knows v + = 0, there is still 
a systematic uncertainty because of the wavelength calibration which was accomplished 
through the OH Meinel band. From their stated "statistical scatter" of 0.5 km s _1 we judge 
that the systematic error in -^ 2 d 3/2 - 2 p 3/2 is no larger than 0.023 cm" 1 (see equation 2), and 
0.012 A in the wavelengths of the red lines. Based on the model-fitting uncertainties alone, 
our value 13637.34 ± 0.01 ±0.03 might seem significantly different, but the values are in fact 
consistent when one accounts for the systematic uncertainties. If we adopted the Sharpee 
et al. (2004) value of E2 Da/2 _2p 3/2 , the velocities in Table 5 would be more consistent for 
37W, but correspondingly less consistent for the other two positions. 

To complete their set of O II 2p 3 energy levels, the blue lines from three nebulae (Bowen 
1960; Baldwin et al. 2000; Sharpee et al. 2004) were used to deduce £ 4S3/2 _2 D5/2 = 26810.76 ± 
0.08 (no systematic error was estimated). This energy difference compares well with our 
result using 23 planetary nebulae, 26810.77 ± 0.03 ± 0.03. 

The UV wavelengths depend on both -E 4 s 3/2 - 2 d V2 and -E' 2 d 3/2 - 2 p 3/2 ) and so the Sharpee 
et al. (2004) values A2470.343± 0.005 and A2470.220± 0.005 are consistent with ours though 
slightly less accurate. 
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6. Conclusions 

From a detailed analysis of Orion nebula and planetary nebulae data, we find new 
energies for the ground configuration. Both model-fitting uncertainties and systematic un- 
certainties are presented. This work confirms the utility of astrophysical measurements in 
determining accurate energies for the O II 2p 3 ground configuration. As a by-product, we 
determine a revised set of air wavelengths of the [O n] visible lines (see Table 5) and the 
UV lines (see § 4). As Sharpee et al. (2004) point out, there are revisions required in the 
standard NIST values. 

Using these revised wavelengths, it is now possible to constrain the velocities of all three 
oxygen zones (O ,O + ,O ++ ), and in turn, the spatial and physical origin of oxygen permitted 
lines can be inferred (Blagrave & Martin 2004b). 

This work was supported by the Natural Sciences and Engineering Research Council 
of Canada. Line wavelengths were obtained from the Atomic Line List 1 maintained by 
P. A. M. van Hoof. The authors wish to thank J. A. Baldwin for obtaining the spectra, and 
G. J. Ferland and E. M. Verner for their comments on this paper. 
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Fig. 1. — Grotrian diagram of II ground configuration energy level transitions (not to 
scale). Wavelengths are approximate air values. 
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Fig. 2. — Gaussian profile fit of all six visible II lines (1SW). The dark solid line repre- 
sents the observed line profile and the lighter dashed line represents the Gaussian fit. The 
wavelength has been adjusted to the rest-frame velocity of the H + gas. Refer to Table 3 for 
reddening-corrected fluxes and the results of the line-fitting. 
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Fig. 3. — Gas velocity as a function of the emitting species' ionization potential. The three 
lines-of-sight have been offset in the x-direction to properly show the standard deviations 
about the mean velocities for each line-of-sight. The ionization potentials of O and ++ 
are 13.6181 eV and 54.934 eV, respectively. + ionization potential is 35 eV. 



Table 1. O II ground configuration energy levels (cm 1 ) 



Designation 


Bowen 


De Robertis 


Level 


Eriksson 


Martin et al. 


This Work a 


Level 




(1955) 


et al. (1985) 


Difference 


(1987) 


(1993) 




Difference a 


4 S3/2 

2 D 5 /2 
2 D 3 /2 

2 P 3 /2 
2 Pl/2 


0.0 
26810.7 
26830.5 

40468.3 


0.0 
26810.5 
26830.6 
40468.1 
40470.1 


26810.5 
20.1 ±0.1 
13637.5 
2.00 ±0.03 


0.00 
26810.52 
26830.57 
40467.69 
40469.69 


0.00 
26810.55 
26830.57 
40468.01 
40470.00 


0.00 

26810.77 ± 0.03 ± 0.03 
26830.57 ±0.01 ±0.03 
40467.91 ±0.02 ±0.05 
40469.93 ± 0.02 ± 0.05 


26810.77 ± 0.03 ± 0.03 
19.80 ±0.01 ±0.00 
13637.34 ±0.01 ±0.03 
2.02 ±0.01 ±0.01 



uncertainties are model-fitting and systematic, respectively 



Table 2. Determination of 2 D 5 / 2 energy from planetary nebulae data 



Nebula 4 S 3 /2- 2 D 3/2 (A) 4 S 3/2 - 2 D 5/2 (A) Weight Reference 



NGC 2440 


3726.00 


3728.69 


0.25 


Hyung & Aller (1998) 


NGC 6543 


3726.06 


3728.81 


1.0 


Hyung ct al. (2000) 


NGC 6567 


3726.25 


3728.97 


0.5 


Hyung ct al. (1993) 


NGC 6572 


3726.01 


3728.76 


1.0 


Hyung ct al. (1994b) 


NGC 6741 


3726.15 


3728.90 


1.0 


Hyung & Aller (1997a) 


NGC 6790 


3726.00 


3728.74 


1.0 


Aller ct al. (1996) 


NGC 6818 


3726.32 


3728.79 


0.0 


Hyung ct al. (1999b) 


NGC 6884 


3726.02 


3728.80 


0.5 


Hyung ct al. (1997) 


NGC 6886 


3726.05 


3728.78 


1.0 


Hyung ct al. (1995) 


NGC 7009major 


3725.97 


3728.71 


1.0 


Hyung & Aller (1995b) 


NGC 7009minor 


3726.14 


3728.95 


0.25 


Hyung & Aller (1995a) 


NGC 7662 


3725.94 


3728.69 


1.0 


Hyung & Aller (1997b) 


IC 351 


3726.04 


3728.81 


1.0 


Feibelman ct al. (1996) 


IC 418 


3726.04 


3728.80 


1.0 


Hyung ct al. (1994a) 


IC 2149 


3725.99 


3728.76 


1.0 


Feibelman ct al. (1994) 


IC 2165 


3726.18 


3728.96 


0.5 


Hyung (1994) 


IC 4634 


3726.27 


3729.02 


1.0 


Hyung ct al. (1999a) 


IC 4846 


3726.17 


3728.93 


1.0 


Hyung ct al. (2001c) 


IC 4997 


3725.97 


3728.71 


1.0 


Hyung ct al. (1994c) 


IC 5117 


3726.03 


3728.75 


0.5 


Hyung ct al. (2001a) 


IC 5217 


3726.04 


3728.77 


1.0 


Hyung ct al. (2001b) 


BD +30 3639 


3726.01 


3728.77 


1.0 


Aller & Hyung (1995) 


Hubble 12 


3726.04 


3728.80 


1.0 


Hyung & Aller (1996) 


Hu 1-2 


3726.07 


3728.83 


1.0 


Pottasch ct al. (2003) 
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Table 3. Best-fit line parameters for visible wavelength transitions. 



Transition 


Position 


Measured Wavelength (A) 


FWHM (km s- 1 ) 


Reddening corrected line strength 






(10~ 13 erg s _1 cm~ 2 ) 


4 S3/2 _ 2 D3/2 


1SW 


3726.072 ± 0.001 


16.8 ±0.2 


426 ±6 


x2 


3726.078 ± 0.001 


18.0 ±0.3 


273 ±4 




37W 


3726.060 ± 0.002 


15.3 ±0.3 


466 ±9 


4 S3/2 _ 2 Ds/2 


1SW 


3728.824 ± 0.001 


17.3 ±0.2 


188 ±2 


x2 


3728.829 ± 0.002 


19.1 ±0.2 


138 ±2 




37W 


3728.811 ± 0.002 


15.7 ±0.3 


213 ±4 


2 D 5 /2- 2 Pl/2 


1SW 


7319.173 ± 0.006 


18.0 ±0.6 


14.7 ±0.5 


x2 


7319.178 ± 0.008 


21.3 ±0.8 


6.9 ±0.2 




37W 


7319.099 ± 0.011 


16.9 ± 1.1 


10.8 ±0.6 


2 D 5/2 - 2 P 3/2 


1SW 


7320.253 ± 0.002 


17.5 ±0.2 


42.6 ±0.5 


x2 


7320.250 ± 0.003 


19.3 ±0.3 


17.1 ±0.2 




37W 


7320.181 ± 0.004 


16.9 ±0.4 


32.0 ±0.6 


2 D3/2- 2 Pl/2 


1SW 


7329.787 ± 0.004 


17.9 ±0.4 


21.1 ±0.4 


x2 


7329.798 ± 0.006 


23.4 ±0.6 


8.6 ±0.2 




37W 


7329.725 ± 0.005 


16.0 ±0.5 


16.6 ±0.4 


2 D3/2~ 2 P3/2 


1SW 


7330.886 ± 0.003 


18.1 ±0.3 


24.0 ±0.4 


x2 


7330.885 ± 0.004 


19.8 ±0.4 


9.5 ±0.2 




37W 


7330.818 ± 0.005 


17.7 ±0.5 


18.2 ±0.5 



Table 4. Comparison between observed and predicted line ratios as determined from transition probabilities, A^. 



Transition 


1SW 




x2 


37W 


Zeippen (1987) 


Wiese et al. (1996) 




Flux 


Ratio 


Flux 


Ratio 


Flux 


Ratio 


Aij 


Ratio 




Ratio 


2 D 5 /2- 2 Pi/2 (A7319) 
2 D 3 /2- 2 Pi/2 (A7330) 


14.7 ±0.5 
21.1 ±0.4 


0.70 ±0.03 


6.9 ±0.2 
8.6 ±0.2 


0.80 ± 0.03 


10.8 ±0.6 
16.6 ±0.4 


0.65 ± 0.04 


0.0563 
0.0941 


0.598 


0.0519 ±0.0052 
0.0867 ± 0.0087 


0.599 ± 0.085 


2 Ds/2- 2 P3/2 (A7320) 
2 D 3 /2- 2 P 3 /2 (A7331) 


42.6 ±0.5 
24.0 ±0.4 


1.77 ±0.04 


17.1 ±0.2 
9.5 ±0.2 


1.80 ±0.04 


32.0 ±0.6 
18.2 ±0.5 


1.76 ±0.06 


0.1067 
0.0580 


1.84 


0.0991 ± 0.0099 
0.0534 ± 0.0053 


1.86 ±0.26 
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Table 5. Line velocities in Orion using recent wavelength tabulations. 



Air A (A) Velocity d (km s _1 ) 



Transition 


Martin a 


Eriksson b 


This work c 


Position 


Observed A (A) 


Martin 


Eriksson 


This work 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


4 S3/2 — 2 Da/2 


3726.032 


3726.032 


3726.032 


1SW 


3726.072 


3.2 ±0.1 


3.2 ±0.1 


3.2 ± 0.1 






±0.001 


x2 


3726.078 


3.7 ±0.1 


3.7± 0.1 


3.7 ±0.1 








±0.004 


37W 


3726.060 


2.3 ±0.2 


2.3 ±0.2 


2.3 ±0.2 


4 S3/2 _ 2 D5/2 


3728.815 


3728.819 


3728.784 


1SW 


3728.824 


0.7 ±0.1 


0.4 ±0.1 


3.2 ± 0.1 






±0.004 


x2 


3728.829 


1.1 ±0.2 


0.8 ±0.2 


3.6 ± 0.2 








±0.004 


37W 


3728.811 


-0.3 ± 0.2 


-0.6 ± 0.2 


2.2 ± 0.2 


2 Ds/2- 2 Pl/2 


7318.92 


7319.07 


7319.073 


1SW 


7319.173 


10.4 ±0.2 


4.2 ±0.2 


4.1 ± 0.2 






±0.009 


x2 


7319.178 


10.6 ±0.3 


4.4 ±0.3 


4.3 ± 0.3 








±0.012 


37W 


7319.099 


7.3 ±0.5 


1.2 ±0.5 


1.1 ± 0.5 


2 D5/2- 2 P 3 /2 


7319.99 


7320.14 


7320.157 


1SW 


7320.253 


10.8 ±0.1 


4.6 ±0.1 


3.9 ±0.1 






±0.009 


x2 


7320.250 


10.6 ±0.1 


4.5 ±0.1 


3.8 ± 0.1 








±0.012 


37W 


7320.181 


7.8 ±0.2 


1.7 ± 0.2 


1.0 ±0.2 


2 D 3 /2- 2 Pl/2 


7329.67 


7329.83 


7329.699 


1SW 


7329.787 


4.8 ±0.1 


-1.8 ±0.1 


3.6 ±0.1 






±0.005 


x2 


7329.798 


5.2 ±0.2 


-1.3 ±0.2 


4.0 ±0.2 








±0.012 


37W 


7329.725 


2.2 ±0.2 


-4.3 ±0.2 


1.1 ± 0.2 


2 D3/2— 2 ?3/2 


7330.73 


7330.91 


7330.786 


1SW 


7330.886 


6.4 ±0.1 


-1.0 ±0.1 


4.1 ± 0.1 






±0.005 


x2 


7330.885 


6.3 ±0.2 


-1.0 ±0.2 


4.0 ±0.2 








±0.012 


37W 


7330.818 


3.6 ±0.2 


-3.8 ±0.2 


1.3 ± 0.2 



a Martin et al. (1993) 
b Eriksson (1987) 

c uncertainties are model-fitting and systematic, respectively, from energy level differences (see Table 1) 
Uncertainties are simply from the least-squares Gaussian fit for each line. 



